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Detailed comparisons are reported between laboratory observations of electron- 
scale dissipation layers near a reconnecting X-line and direct two-dimensional 
full-particle simulations. Many experimental features of the electron layers, 
such as insensitivity to the ion mass, are reproduced by the simulations; the 
layer thickness, however, is about 3 — 5 times larger than the predictions. 
Consequently, the leading candidate 2D mechanism based on collisionless elec- 
tron nongyrotropic pressure is insufficient to explain the observed reconnec- 
tion rates. These results suggest that, in addition to the residual collisions, 
3D effects play an important role in electron-scale dissipation during fast re- 
connection. 
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Despite the disruptive influences of magnetic reconnection on large-scale structures in 
plasmas, the crucial topological changes and their associated dissipation take place only 
within thin current layers. The classical collisional models, where electrons and ions flow 
together through a single thin and long layer, fail to explain the observed fast reconnection 
rates. Modern collisionless models predict [Sonnerup, 1979; Mandt et al, 1994; Birn 
et al, 2001] that ions exhaust through a thick, ion-scale layer while mobile electrons 
flow through a thin, electron-scale layer, allowing for efficient release of magnetic energy. 
These ion layers have been frequently detected in space [e.g. Deng and Matsumoto, 2001; 
0ieroset et al, 2001; Mozer et al, 2002] and studied in detail in the laboratory [Ren 
et al, 2005; Yamada et al, 2006; Brown et al, 2006]. In contrast, the electron layers, 
where magnetic field dissipates, are rarely encountered in space and are often detected at 
places far from the reconnection X-line line [Scudder et al, 2002; Mozer, 2005; Wygant 
et al, 2005; Phan et al, 2007]. Therefore, whether the electron layers indeed exist near 
the X-line, and if yes, whether their associated dissipation results predominantly from 
laminar two-dimensional (2D) or three-dimensional (3D) dynamics as suggested by Xiao 
et al. [2006, 2007], is still an open question. Here we report detailed comparisons between 
recent laboratory observations of the electron layers near the X-line [Ren et al, 2008] 
and direct full-particle simulations in 2D. The measured electron layers display properties 
strikingly similar to predictions by 2D particle simulations, including their geometrical 
shape, insensitivity to ion mass, and sensitivity to the boundary conditions, but disagree 
on the electron layer thickness. As a consequence, the leading 2D mechanism based on 
collisionless electron nongyrotropic pressure is shown to be largely insufficient to explain 
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the observed reconnection rates. These results suggest that, in addition to the residual 
Coulomb collisions, 3D effects play an important role in electron-scale dissipation during 
fast reconnection. 

The laboratory measurements were performed on the well controlled and diagnosed 
experiment, Magnetic Reconnection Experiment (MRX) [Yamada et al, 1997], as il- 
lustrated in Fig.l. A pair of coil assemblies, known as flux-cores, are used to axisym- 
metrically initiate and maintain the reconnection process. Plasma is made by ioniz- 
ing a pre-filled gas through pulsing toroidal field coil current within the flux-cores dur- 
ing the period when the current flowing in the poloidal field (PF) coils peaks. When 
the PF coil current is ramped down after the plasma is made, the field lines wrapped 
around both flux-cores are "pulled" back, reconnect, and move towards the flux-cores. 
Most of the important quantities can be either directly determined or indirectly inferred 
from these measurements in cylinderical coordinates (R, Z, 6) assuming axisymmetry: 
poloidal flux ip(R,Z,t) = J^27iR'Bz{R',Z,t)dR' where Bz is the reconnecting field; the 
toroidal reconnection electric field E e = (dip/dt)/2nR; and the toroidal current density 
je ~ n^dBz/dR. The density n and electron temperature T e are measured by a triple 
Langmuir probe and the flow speeds are determined by a Mach probe. The typical plasma 
parameters are: n ~ (0.1 - 2) x 10 20 m" 3 , T e ~ T { ~ (3 - 15) eV, B < 0.5 kG. 

Detection of the electron dissipation layer is made possible by taking advantage of the 
differential motions between electrons and ions or the so-called Hall effects [Sonnerup, 
1979] in the reconnection region without a guide field. These differential motions (or elec- 
tric current) within the reconnection plane produce out-of-plane magnetic field component 
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{Be) with a quadrupole shape. Conversely, accurate measurements of the Be profile can 
determine the in-plane electron flow because of the much slower ion flow in this region, 
and thus characterize the electron dissipation layer. These measurements are performed 
using five linear arrays of pickup coils (Fig.l); each array measures a one- dimensional 
profile of Be with a frequency response of 300kHz and with spatial resolutions up to 2.5 
mm. This distance is close to the electron skin depth, c/u pe (=0.7-1. 5mm) where u pe is 
the electron plasma angular frequency, and adequately resolves the electron layer whose 
minimum full thickness is 10 mm (see below). These arrays are housed by thin glass 
tubes of outer diameter of 4 mm (four arrays) or 5 mm (one array) with shielding from 
electrostatic noise. The presence of these probes in the plasma does not appear to affect 
the reconnection process, but it may cause modest overestimates of the electron layer 
thickness (see below). 

One such example measurement is shown in Fig. 2(b) where the in-plane electron flow 
(V e z and V c r) is shown as arrows while the normalized, out-of-plane magnetic field is 
shown as color-coded contours in the left half of the reconnection plane. Electron outflow 
speed, V e z, is also shown as functions of Z in Fig. 2(c) (at the current sheet center) and 
R in Fig. 2 (a) (across the reconnection region at the location where V eZ peaks). The 
dimensions of the electron layer can be characterized by the half thickness 5 e (the radial 
distance during which V e z decreases by 60% from its peak value) and the half length 
L e (the axial distance during which V eZ increases from zero to its peak). We positively 
identify this region as the electron dissipation layer because both its dimensions, 5 e and 
L e , are independent of ion mass, as shown in Fig. 3 for 5 e . 
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Dissipation in the electron layer is governed by the electron equation of motion, 

m e n (J^ + V e- V ) V e = -en (E + V e x B) - V ■ P e + enrjspitzerj , (1) 

where m e is electron mass, P e electron pressure tensor, and r\ Spitzer the Spitzer resistivity 
due to Coulomb collisions with ions [Spitzer, 1962]. In the modern collisionless steady- 
state 2D models, the reconnection electric field, Eg, can be only possibly balanced by 
either the Hall term (V e x B) e (j x B)g/en, the inertia terms, or the electron pressure 
tensor term (V • P e )e- While the Hall term is important in supporting Eg within the ion 
layer [Birn et al, 2001; Ren et al, 2005; Yamada et al, 2006], it diminishes within the 
electron layer especially near the X-line due to the vanishing B. It has been shown in 
particle simulations [Cai and Lee, 1997; Hesse et al, 1999; Pritchett, 2001; Kuznetsova 
et al, 2001] that Eg is supported primarily by the electron pressure tensor term following 
earlier suggestions [Vasyliunas, 1975]. This mechanism has been since widely accepted as 
the leading candidate to provide the required dissipation within the electron layer. It is, 
however, extremely difficult to confirm this pressure anisotropy, directly or indirectly, by 
measurements in real plasmas [Scudder et al, 2002]. 

One of the predictions of these 2D particle simulations is that the half thickness of 
the electron layer, S e , scales as (1 — 2)c/u pe [Pritchett, 2001]. The measured S e in MRX, 
however, scales as ~ 8c/w pe (Fig.3). Current blockage due to the probes is estimated 
to lead to a 6 — 44% increase in the measured 5 e , depending on the ratio of 5 e to the 
glass tube radius. Applying these corrections leads to 5 e = (5.5 — 7.5)c/u pe . To better 
compare with the experiment, on the other hand, we have constructed a kinetic numerical 
model [Dorfman et al, 2008] using boundary conditions similar to the MRX based on 
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the existing NPIC 2D code [Daughton et al, 2006]. A 75cmx 150cm simulation box is 
used with conducting boundary conditions for fields and elastic reflection for particles 
at the walls. Two current carrying coils of radius 1.3 cm are contained within a larger 
concentric flux core of radius 9.4 cm. The flux cores are spaced 40 cm apart as in the 
experiment. The flux core surface is approximated as an insulating boundary; particles 
may be absorbed or reflected. Due to constraints on computation resources, the number 
of the Debye lengths per c/u pe is limited compared to the experiment, but there is strong 
evidence that the reconnection rate and electron layer scalings are insensitive to this 
number as long as the initial plasma beta is fixed [Dorfman et al, 2008]. As the current 
is ramped down according to a sinusoidal waveform modeled on the PF coil current of 
MRX and reconnection is driven, both ion and electron dissipation layers are formed. 
Simulation parameters are chosen such that the global reconnection rate and the current 
sheet thickness on the ion scale match the observations. An example run is shown in 
Fig.2(d-f) in the same format as in Fig.2(a-c), and most of the observed features, including 
geometrical shapes and out-of-plane magnetic component, are reproduced. 

The quantitative agreement between experiment and simulation is, however, found for 
only the global ion dynamics but not the local electron dynamics. The experimentally 
observed independence of 5 e and L e on ion mass was reproduced as shown for 5 e by the 
open squares in Fig. 3 for a fixed but artificially heavy electron mass. The values of 5 e 
in units of c/w pe (evaluated using a line-averaged density at Z — 0), however, are much 
smaller in simulations than in experiments, as illustrated by an alternative ordinate in 
Fig.2(a) and (d). In Fig. 3, a case at higher mass ratio (400) with a different electron mass 
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is also plotted along with simulations with a realistic hydrogen mass ratio but a smaller 
simulation domain and open boundary conditions [Daughton et ai, 2006]. All of these 
cases, including more recent simulations using different open boundary conditions [Huang 
and Ma, 2008], confirm a linear relation of 5 e = (f .5 — 2)c/u pe which is about 3 — 5 times 
thinner than the experiment. In contrast, the dependence of the length of the electron 
layer (L e ) on c/u pe is less robust; it can change significantly when the reflection coefficient 
parameter on the flux core surface is varied [Dorfman et ai, 2008] as expected from the 
observed dependence of the reconnection process on boundary conditions [Kuritsyn et ai, 
2007]. 

The fact that the observed electron layers are substantially thicker than the numerical 
predictions implies different dissipation mechanisms operating between these two cases. 
In fact, our collisionless simulation model does not include the residual collisions between 
electrons and ions or neutrals. But in MRX only a fraction of Eg can be accounted 
for by the classical resistivity, E v = r)s P itzerje (Fig.4). Collisions between electrons and 
neutrals, and electron collisional viscous effects are also estimated to be unimportant in 
these discharges with low fill pressure. The electron inertia terms, (m e /e)[(V e Rd/dR) + 
(V eZ d/dZ)]V ee , are estimated to be on the order of 1 V/m, which is negligibly small. Near 
the X-line, the effects due to electron nongyrotropic pressure can be well approximated 
by [Hesse et ai, 1999] 




as also validated in our kinetic model. Direct evaluations of E NG using the measured pro- 
file, V e z(Z) as in Fig. 2, gives values only a small fraction of E e — E v (Fig. 4). This leaves 
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the majority of E$ still unexplained, and therefore there must exist additional dominant 
dissipation mechanisms. Because our kinetic model contains all possible collisionless ki- 
netic mechanisms operative in 2D, these dominant mechanisms must be 3D in character, 
including effects due to current sheet deformation or plasma turbulence through wave- 
particle interactions within the current sheet. The latter was indeed already suggested by 
the detection of electromagnetic fluctuations [Ji et al, 2004] when dissipation increases 
at low collisionalities [Ji et al, 1998]. This subject is also under intensive theoretical 
and numerical investigation, such as recently by Moritaka et al. [2007], in the search for 
mechanisms for fast reconnection. Lastly, we comment that these 3D effects, in additional 
to the residual collisions, may diffuse substantially the predicted two-scale structures seen 
in the profiles of the reconnecting magnetic field, which remain undetected thus far in the 
experiment. 
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Flux cores 



Figure 1. Experimental set-up of MRX device. The toroidal direction points along the 
current sheet while the poloidal direction wraps around the flux cores. 
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Figure 2. Identification of electron dissipation layer. The top three panels (a-c) show 
an experimental example taken from a hydrogen plasma with a fill pressure of 2 mTorr. 
Results from a corresponding simulation are shown in the same format in the lower three 
panels (d-f). The parameters used in the simulation are: 864 x 1728 cells with 0.5 billion 
particles per species, initial density of 2.6 x 10 19 m~ 3 , m; = m hydrogen , m e = m hydrogen /75, a 
time scale for the coil current ramp down is 185 initial ion cyclotron times, and no particle 
reflections at the flux core surface. 
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Figure 3. Scaling of width of electron dissipation layer. Filled symbols show the 
experimentally measured 5 e as a function of the electron skin depth (c/u pe ) for three 
different ion species. The error bars result mainly from shot-to-shot variations. Open 
symbols show S e determined from 2D PIC simulations. 
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Figure 4. Composition of reconnecting electron field, Eg, for (a) helium and (b) 
deuterium plasmas. Total reconnecting electric field in MRX, Eg, and the part of it due 
to electron-ion collisions, E v = rjspUzerje near the X-line are plotted as a function of c/u pe . 
The estimated electric field due to electron nongyrotropic pressure, E NG , is also shown. 
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